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shears  of  the  zonal  wind.  The  individual  cyclones  which  made  up  the 
composite  study  were  examined  for  these  same  parameters  and  the  results 
obtained  were  also  significant. 

Based  on  the  results  obtained  from  testing  the  individual  storm 
systems,  a  scheme  was  developed  for  forecasting  intensity  change  of 
tropical  cyclones.  The  forecasting  scheme  was  tested  against  independ¬ 
ent  storm  systems  and  was  found  to  give  correct  predictions  of  inten¬ 
sification  or  non- intensification  approximately  75%  of  the  time.  The 
greatest  strength  of  the  scheme  is  that  it  gives  a  quantitative  index 
which  can  be  obtained  from  purely  objective  techniques. 


The  technique  was  also  applied  to  other  comparative  rawinsonde 
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able  to  improve  our  skill  at  forecasting  future  tropical  cyclone 
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ABSTRACT  OF  THESIS 


TROPICAL  CYCLONE  INTENSITY  CHANGE  -  A  QUANTITATIVE 
FORECASTING  SCHEME 

This  paper  discusses  tropical  cyclone  intensity  change  based  on  the 
results  obtained  from  an  observational  study  comparing  tropical  systems 
which  intensify  and  those  that  do  not.  Tropical  cyclones  found  in  the 
Gulf  of  Mexico  during  the  period  1957-1977  were  studied  because  of  the 
abundance  of  data  in  this  region  and  also  for  investigation  of  the 
effects  of  large  scale  baroclinic  interaction  on  the  northwest  side  of 
the  cyclone. 

Initially  weather  charts  were  examined  in  hopes  of  finding  comnon 
meteorological  parameters  which  were  different  between  the  intensifying 
and  non- intensifying  storm  systems.  But  when  obvious  differences  could 
not  be  found  from  the  charts,  a  rawinsonde  composite  analysis  was  per¬ 
formed  to  investigate  mean  differences  between  the  two  systems.  The 
results  obtained  wer significant.  Very  large  differences  were  found  on 
the  northern  side  of  the  cyclone  particularly  in  the  200  mb  height  field, 
the  850  mb  temperature  field,  the  200  mb  zonal  wind  and  the  vertical 
shears  of  the  zonal  wind.  The  individual  cyclones  which  make  up  the 
composite  study  were  examined  for  these  same  parameters  and  the  results 
obtained  were  also  significant. 

Based  on  the  results  obtained  from  testing  the  individual  storm 
systems,  a  scheme  was  developed  for  forecasting  intensity  change  of 
tropical  cyclones.  The  forecasting  scheme  was  tested  against  independ¬ 
ent  storm  systems  and  was  found  to  give  correct  predictions  of  inten¬ 
sification  or  non- intensification  approximately  75%  of  the  time.  The 
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greatest  strength  of  the  scheme  is  that  it  gives  a  quantitative  index 
which  can  be  obtained  from  purely  objective  techniques. 

The  technique  was  also  applied  to  other  comparative  rawinsonde 
composite  studies  done  in  other  ocean  basins.  Again,  the  results  were 
highly  significant.  There  appears  to  be  a  definite  relationship  between 
intruding  baroclinic  zones  on  tropical  cyclone  systems  and  the  subse¬ 
quent  intensity  change  of  that  system.  By  measuring  those  parameters 
which  are  most  affected  by  baroclinic  interaction,  we  might  be  able  to 
improve  our  skill  at  forecasting  future  tropical  cyclone  intensity 
change. 
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Atmospheric  Science  Department 
Colorado  State  University 
Fort  Collins,  Colorado  80523 
Spring,  1981 
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1.  INTRODUCTION 


Forecasting  the  intensity  change  of  tropical  cyclones  has  eluded 
forecasters  for  years.  Despite  recent  advances  in  computer  and  satellite 
technologies  the  forecasting  skill  has  not  shown  much  improvement 
(Hebert,  1978).  The  computer  has  made  more  meteorological  products 
available  and  the  satellite  has  made  us  more  aware  of  the  existence  of 
developing  systems,  but,  to  date  neither  has  led  to  significant  improve¬ 
ments  in  the  general  forecasting  skill. 

In  general,  intensity  change  may  be  due  to  both  internal  and  exter¬ 
nal  mechanisms  such  as  that  given  by  Eq.  1: 


|i  -  A  +  B  (1) 

where  I  =  intensity,  A  =  inner  core  (0-2°  radius)  influences  on  cyclone 
intensity  change,  and  B  =  outer  region  (2-10°  radius)  influences  on 
cyclone  intensity  change.  Limited  observations  of  the  inner  core  do 
not  allow  us  to  accurately  measure  the  physical  parameters  which  may 
lead  to  intensification.  However,  new  evidence  (Nunez,  1981)  indicates 
that  the  large  scale  environment  may  have  a  significant  effect  on  the 
intensification  process.  Therefore,  this  study  examines  the  surround¬ 
ing  outer  region  environmental  influences  (effect  B  in  Eq.  1)  to 
determine  how  they  bring  about  cyclone  intensity  change. 

The  research  findings  of  Nunez  (1981)  on  tropical  cyclone  intensity 
change  formed  a  background  for  this  study.  Nunez  examined  cyclone  in¬ 
tensity  change  by  comparing  composite  cases  of  deepening  versus  filling 
tropical  cyclones  in  the  northwest  Pacific  and  northwest  Atlantic.  He 
found  that  filling  systems  had  stronger  upper  level  height  and 
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temperature  gradients  on  the  north  side  of  the  storm,  and  often  an  upper 
level  trough  to  the  northwest  of  the  storm  center.  He  attributed  the 
differences  between  the  filling/deepening  data  sets  to  strong  baroclinic 
interaction  on  the  north  and  northwest  octants  of  the  storm.  His  hypo¬ 
thesis  was  that  an  imbalance  exists  in  the  thermal  wind  equation  such 
that  the  deepening  systems  have  a  larger  vertical  shear  of  the  tangential 
wind  (WS)  than  would  be  expected  from  the  baroclinicity  (B) .  As  a 
result  of  this  imbalance,  the  inner  area  of  the  disturbance  warms  through 
vertical  mass  recycling  which  produces  intensification.  The  imbalance 
continues  as  long  as  cumulus  clouds  and  vertical  mass  recycling  occur 
in  the  favorable  vertical  shear  areas.  Under  unfavorable  vertical  shears, 
cumulus  activity  acts  to  reduce  vertical  wind  shear  causing  a  decrease 
in  the  baroclinicity  and  weakening  of  the  cyclone. 

Based  on  the  Nunez  composite  data  results,  the  author  has  under¬ 
taken  the  investigation  of  such  baroclinic  interaction  on  an  individual 
case  basis.  Tropical  cyclones  in  the  Gulf  of  Mexico  were  selected  for 
investigation  because  of  the  abundance  of  data  on  the  north  side  of  the 
storm.  The  storms  were  stratified  as  intensifying  and  non- intensifying. 

It  was  hoped  that  the  effects  of  large  scale  baroclinic  interaction  on 
the  north  side  of  the  non-intensifying  cyclone  systems  would  produce 
evidence  of  intensity  change  that  could  be  recognized  from  the  weather 
charts.  But  after  hours  of  detailed  analysis  of  upper  and  lower  tropo¬ 
spheric  weather  charts,  the  author  was  convinced  that  the  charts  did  not 
provide  obvious  observational  differences  which  would  signal  an  intensity 
change.  The  presence  of  trough  interaction  was  observed  but  it  was  dif¬ 
ficult  to  quantify  the  effects  of  that  interaction  solely  from  compari¬ 
sons  of  the  charts.  Therefore,  rawinsonde  composites  were  formed  for 


both  data  sets  and  the  techniques  developed  by  Gray  (1967) ,  Williams  and 
Gray  (1973),  Ruprecht  and  Gray  (1976),  Zehr  (1976),  Erickson  (1977), 
Arnold  (1977),  McBride  (1979),  and  Nunez  (1981)  were  used  to  determine 
the  major  physical  differences  between  systems.  The  knowledge  gained 
from  this  comparative  composite  study  was  then  applied  to  individual 
case  data  to  develop  a  purely  objective  and  quantitative  scheme  for 
forecasting  tropical  cyclone  intensity  change.  Independent  data  were 
selected  to  further  test  the  forecast  scheme  and  finally  a  thorough 
statistical  analysis  was  made  of  all  data  sampled. 

1.1  Compositing  Philosophy 

The  Gulf  of  Mexico  is  one  of  the  best  locations  in  the  world  to 
study  individual  storms  because  of  the  abundance  of  surface  and  upper 
air  observing  sites  throughout  the  Caribbean  and  southeastern  United 
States  (Fig.  1).  But,  even  in  the  Gulf  of  Mexico  problems  arise  be¬ 
cause  of  the  sparcity  of  data  over  Mexico  and  the  lack  of  observations 
within  the  center  of  the  Gulf.  The  routine  observational  data  is  often 
augmented  by  aircraft  observations  but  these  data  are  often  limited  to 
inner  core  information  and  rarely  at  more  than  one  or  two 
levels.  The  compositing  technique  allows  us  to  overlook  these  short¬ 
comings  by  using  rawinsonde  data  for  many  cyclones  which  show  similar 
characteristics. 

The  composite  technique  tends  to  smooth  over  many  of  the  individual 
characteristics  and  diurnal  variations  of  a  single  storm.  Still,  by 
proper  data  selection  and  handling  techniques,  we  can  gain  valuable 
information  despite  the  loss  of  some  of  these  'eddy  qualities'. 

This  study  takes  the  general  rawinsonde  composite  approach  to  de¬ 
termine  differences  between  the  two  data  sets  described  below.  The 


Fig.  1.  West  Atlantic  rawinsonde  data  network  utilized  in  the  present 
study. 

results  are  then  applied  to  the  individual  cases  to  see  how  representa¬ 
tive  the  composite  differences  are. 


1.2  Compositing  Technique 

The  composite  is  accomplished  by  positioning  a  cylindrical  grid 
with  a  radius  of  15°  latitude  at  the  surface  center  of  the  storm  for 


each  time  period.  The  grid  consists  of  21  vertical  levels  extending 
from  the  surface  to  50  mb.  The  horizontal  grid  is  divided  into  eight 
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octants  and  eight  radial  belts  as  depicted  in  Fig.  2.  This  subdivision 
yields  64  boxes  whose  areal  extent  increases  radially  outward. 

By  using  this  grid  system  over  every  storm  in  our  set,  we  can 
accumulate  many  values  of  each  parameter  in  each  of  the  64  boxes.  The 
accumulated  values  in  each  grid  box  are  averaged  and  the  value  is 
assigned  to  represent  the  corresponding  box.  Repeating  this  averaging 
process  for  each  of  the  meteorological  parameters  in  each  of  the  64  grid 
boxes  yields  the  basic  structure  for  the  composited  cyclone. 

The  observed  or  computed  parameters  include  14  dynamic  and  10 
thermodynamic  parameters  as  listed  below-: 


Dynamic  Parameters 

u  (zonal  wind) 

v  (meridional  wind) 

(radial  wind) 

V  (tangential  wind) 

u 

Vp  (wind  parallel  to  storm  motion) 

(wind  normal  to  storm  motion) 

Vrff  (Coriolis  torque) 

(horizontal  moisture  transport) 

DIV  (divergence) 

w  (vertical  motion) 

q  m  (vertical  moisture  transport) 
c 

c;  (vortlcity) 

3u  (vertical  zonal  wind  shear) 

3p 


Thermodynamic  Parameters 

T  (temperature) 

Tv  (virtual  temperature) 

0  (potential  temperature) 

Z  (height) 

RH  (relative  humidity) 

q  (specific  humidity) 

s  (static  energy) 

h  (moist  static  energy) 

h*  (saturated  moist  static  energy) 

3gz  (radial  height  gradient) 

3r 


All  parameters  were  examined  in  this  study  but  only  the  zonal  wind,  verti 
cal  shear  of  the  zonal  wind,  height  and  temperature  were  found  to  exhibit 
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Fig.  2.  Compositing  grid.  Arrow  points  north. 

strong  differences  between  the  intensifying  and  non- intensifying  cyclone 
systems. 

Corrections  in  the  relative  position  of  the  rawinsonde  to  the  storm 
center  are  made  by  assuming  the  balloon  is  released  30  minutes  prior  to 
the  scheduled  observation  time  and  that  the  ascent  rate  is  5  m  s  . 

Drift  corrections  are  then  made  by  using  the  reported  wind  profiles. 
Computational  procedures  also  correct  for  translational  movement  of  the 
storm  while  the  balloon  is  airborne. 

The  computer  programs  are  designed  to  composite  in  four  different 
cylindrical  coordinate  systems,  but  this  study  will  be  limited  to  the 
natural  (NAT)  coordinate  system.  In  this  coordinate  system  the  composi¬ 
ting  grid  is  stationary  with  respect  to  the  storm  center.  Octant  1  is 


always  pointing  to  the  north,  regardless  of  the  storm's  direction  of 


movement. 

1.3  Data  Set 

This  study  utilizes  21  years  (1957-1977)  of  northwest  Atlantic  ra- 
winsonde  data  from  the  stations  shown  in  Fig.  1.  The  data  come  from 
the  Northern  Hemisphere  Data  Tabulations  (NHDT)  tapes  from  Asheville, 
North  Carolina  Records  Center  which  was  made  available  to  the  Colorado 
State  University  (CSU)  project  courtesy  of  Mr.  Roy  Jenne  of  the  National 
Center  for  Atmospheric  Research  (NCAR)  in  Boulder,  Colorado.  Both  00 
and  12  GMT  soundings  were  used. 

The  cyclones  were  selected  from  the  official  best  track  positions 
of  the  National  Hurricane  Center.  The  cyclones  were  chosen  if  they 
made  landfall  along  the  Gulf  Coast  of  the  United  States  from  approxi¬ 
mately  100  miles  south  of  Brownsville,  TX  to  the  Florida  Keys.  Cyclones 
were  stratified  into  intensifying  or  non- intensifying  data  sets  based 
on  their  maximum  wind  speed  changes  for  the  time  they  were  in  the  Gulf. 

A  cyclone  was  considered  an  intensifier  if  it  had  a  wind  speed  increase 
of  greater  than  20  knots  (kts)  in  a  24  hour  period  within  42  hours  of 
landfall  (see  Figs.  3  and  4).  There  were  seven  cyclones  which  were 
classified  as  intensifying  but  showed  a  significant  filling  a  few  hours 
prior  to  landfall.  These  storms  were  not  considered  in  this  study.  A 
complete  listing  of  all  cyclones  which  make  up  the  two  data  sets  is 
shown  in  Table  1.  Table  1  also  lists  the  cyclone's  maximum  sustained 
wind  speed  for  the  LF  -42,  LF  -18  and  landfall  time  periods  (LF  repre¬ 
sents  Land  Fall) .  The  average  intensity  of  the  intensifying  cyclones  is 
nearly  twice  as  strong  as  the  non- intensifying  systems. 


9 


Fig.  4.  Maximum  wind  speed  changes  for  all  storms  which  make  up  the 
non- intensifying  composite  data  set. 

The  composites  were  prepared  for  two  time  periods  -  approximately 
18  and  42  hours  prior  to  the  landfall  period.  These  time  periods  were 
selected  to  approximate  the  critical  warning  periods  used  by  the 
National  Hurricane  Center  (i.e.,  24  and  48  hours  prior  to  landfall). 

The  landfall  period  minus  18  hours  (LF  -18)  was  obtained  by  using  both 
the  LF  -12  and  LF  -24  data  composited  around  the  LF  -18  cyclone  position. 
This  was  done  to  double  the  rawinsonde  reports  because  the  data  samples 
were  small  (22  intensifying  cases  vs.  13  non- intensifying  cases). 
Similarly,  the  LF  -42  composite  was  obtained  by  using  LF  -36  and  LF  -48 
data.  The  number  of  soundings  which  fell  into  each  of  the  64  grid 
boxes  for  the  period  LF  -18  is  shown  in  Fig.  5  for  the  intensifying  and 
non- intensifying  sets.  It  is  easy  to  see  from  this  figure  that  the  most 
abundant  data  sources  lie  to  the  north  of  the  storm  -  a  factor  which  we 
will  take  into  consideration  later.  The  distribution  of  soundings  for 
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TABLE  1 

Listing  of  individual  storms  for  the  intensifying  and  non- intensifying 
composite  sets.  The  maximum  wind  is  indicated  for  the  three  time  periods. 

INTENSIFYING  SET 

MAXIMUM  WIND  (kts) 


STORM 

LF  -42 

LF  -18 

LF 

H1  Audry  25-28  JUN  1957 

60 

80 

125 

T2  Bertha  8-11  AUG  1957 

3 

Missing 

35 

60 

H  Debra  22-27  JUL  1959 

25 

35 

70 

H  Judith  17-21  OCT  1959 

Missing 

30 

65 

H  Carla  3-15  SEP  1961 

110 

130 

150 

H  Cindy  16-19  SEP  1963 

Missing 

30 

70 

T  Abby  5-8  AUG  1964 

25 

25 

55 

H  Isbell  8-16  OCT  1964 

30 

80 

no 

H  Betsy  26  AUG  -  12  SEP  1965 

no 

115 

135 

H  Beulah  5-22  SEP  1967 

95 

110 

140 

H  Abby  1-13  JUN  1968 

30 

60 

65 

T  Candy  22-26  JUN  1968 

Missing 

25 

60 

H  Camille  14-22  AUG  1969 

100 

130 

165 

H  Celia  30  JUL  -  5  AUG  1970 

60 

85 

no 

H  Ella  8-13  SEP  1970 

30 

70 

no 

T  Felice  11-17  SEP  1970 

30 

30 

60 

H  Edith  5-18  SEP  1971 

45 

45 

85 

H  Carmen  29  AUG  -  10  SEP  1974 

70 

105 

130 

H  Caroline  24  AUG  -  1  SEP  1975 

35 

65 

100 

H  Eloise  13-24  SEP  1975 

45 

65 

no 

H  Anita  29  AUG  -  2  SEP  1977 

75 

85 

150 

H  Babe  3-8  SEP  1977 

25 

40 

65 

AVERAGE  INTENSITY 

56 

67 

100 

AVERAGE  CHANGE  OF  MAXIMUM 
SUSTAINED  WIND  SPEED 
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TABLE  1  (cont'd) 
NON -INTENSIFYING  SET 


MAXIMUM 

WIND  (kts) 

STORM 

LF  -42 

LF  -18 

LF 

T  Debbie  7-8  SEP  1957 

Missing 

35 

35 

T  Esther  16-19  SEP  1957 

Missing 

45 

45 

H  Ella  30  AUG  -  6  SEP  1958 

60 

60 

55 

T  Arlene  28  MAY  -  2  JUN  1959 

30 

45 

50 

T  Irene  6-8  OCT  1959 

Missing 

30 

50 

T  TS#1  22-26  JUN  1960 

Missing 

15 

40 

T  Florence  17-26  SEP  1960 

25 

25 

25 

T  TS//1  2-11  JUN  1964 

25 

25 

30 

T  TS//1  11-18  JUN  1965 

25 

35 

45 

H  Gladys  13-21  OCT  1968 

65 

65 

70 

T  Genny  1-6  OCT  1969 

25 

25 

25 

H  Alma  17-27  MAY  1970 

30 

25 

25 

H  Agnes  14-22  JUN  1972 

55 

65 

75 

AVERAGE  INTENSITY 

38 

38 

44 

AVERAGE  CHANGE  OF  MAXIMUM 
SUSTAINED  WIND  SPEED 


TI  represents  a  maximum  intensity  of  hurricane  (>_  65  kts) . 

2 

T  represents  a  maximum  intensity  of  tropical  storm  (>_  35  kts) 

3 

Missing  is  used  because  the  cyclone  had  not  yet  formed. 

4 

TS//1  of  1960  should  have  been  stratified  as  an  intensifying  sys¬ 
tem  based  on  post  analysis.  The  error  could  not  be  corrected  for  the 
composite  study  but  the  system  is  treated  as  an  intensifying  cyclone 
in  sections  3  and  5. 
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I 
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2.  RAWINSONDE  COMPOSITE  COMPARISON:  INTENSIFYING  VS.  NON- INTENSIFY  INC 


Some  comparisons  of  meteorological  parameters  between  the 
intensifying  and  non- intensifying  composite  data  sets  were  made  between 
0-4°  radius  primarily  for  the  investigation  of  structural  differences. 
But  the  majority  of  the  comparisons  were  made  at  outer  radii  where  the 
data  were  most  abundant  and  reliable.  Also,  because  the  intensifying 
data  set  has  an  average  intensity  which  is  nearly  twice  as  large  as  the 
non-intensifying  set,  the  differences  in  parameters  at  close  radii  may 
be  biased  by  this  intensity  difference.  It  is  observed  that  at  outer 
radii  the  differences  in  cyclone  intensity  between  the  two  systems 
becomes  negligible.  Additionally,  the  outer  region  is  where  we  would 
expect  any  baroclinic  interaction  to  be  most  noticeable.  This  section 
examines  the  climatological  differences,  effects  of  movement,  thermo¬ 
dynamic  differences  and  dynamic  differences  for  the  composited  data 
sets. 


2.1  Climatological  Differences 

In  the  Gulf  of  Mexico,  August  5  through  October  20  is  climatologi- 
cally  the  most  active  period  for  tropical  cyclone  development  with  a 
minor  secondary  maximum  occurring  in  June'*'.  Table  2  gives  a  breakdown 
of  the  occurrence  by  month  of  the  cyclones  which  make  up  the  composite 
study.  If  we  let  the  'active  period’  be  the  entire  months  of  August, 
September  and  October  we  see  that  82%  of  the  intensifying  and  54%  of 
the  non- intensifying  cyclones  occur  during  this  period.  If  we  consider 
all  cyclones  from  this  study  which  occur  during  the  active  period  for 

^Determined  from  cumulative  summaries  over  the  years  1886-1977  as 
found  in  Tropical  Cyclones  of  the  North  Atlantic  Ocean,  1871-1977. 
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TABLE  2 

■  r  r  cm  <  .>(  all  storms  during  the  years  1957-1977  which  make 
^ i t •  • i udv  as  listed  in  Table  1. 


MAY  JUN  JUL  AUG  SEP  OCT  TOTAL 


INTENSIFYING 

0 

3 

D 

5 

II 

2 

22 

NON- INTENSIFYING 

2 

□ 

0 

0 

Q 

3 

13 

the  years  1957-1977,  we  see  that  72%  of  them  are  intensifying  systems, 
of  all  cyclones  which  occur  during  May,  June  or  July  only  40%  intensify. 

Preliminary  inspection  of  Table  2  may  lead  to  the  conclusion  that 
the  basic  difference  between  the  intensifying  and  non- intensifying  sys¬ 
tems  results  from  their  date  of  occurrence.  This  is  not  the  case  how¬ 
ever.  Seasonal  climatological  differences  do  occur  but  they  are  not 
dominant.  A  discussion  of  those  differences  will  be  included  in  section 
4.  Allowances  for  these  seasonal  differences  will  be  made  as  the  fore¬ 
casting  scheme  is  developed. 

2.2  Cyclone  Movement 

The  cyclone  tracks  for  both  data  sets  are  shown  in  Figs.  6  and  7. 
These  figures  indicate  that,  in  general,  cyclones  with  a  west  northwest 
(WNW)  track  are  more  likely  to  intensify  while  those  with  a  north  or 
northeast  (NNE)  track  are  more  likely  to  fill.  The  average  direction  of 
movement  from  the  composite  study  is  towards  303°  for  the  intensifying 
data  set  and  towards  330°  for  the  non- intensifying  set  at  the  LF  -42 
time  period.  Averages  for  the  LF  -18  time  period  are  315°  for  the 
intensifying  set  and  342°  for  the  non- intensifying  set.  The  difference 
between  the  data  sets  is  a  consistent  27°  for  both  time  periods. 


NON- INTENSIFYING 
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However,  because  of  erratic  cyclone  movements  in  both  data  sets,  any 
generalizations  about  intensity  change  and  motion  relationships  may  not 
hold  well  for  the  individual  case. 

The  speed  of  movement  between  the  two  data  sets  was  also  analyzed 
and  compared.  In  the  mean  there  was  little  difference.  But  after  a 
cyclone  has  made  a  turn  toward  the  NNE  there  is  a  slight  acceleration  as 
the  cyclone  becomes  caught  up  in  the  westerly  wind  zone.  A  correlation 
could  not  be  found  between  the  rate  of  deepening  or  filling  and  the 
speed  of  movement. 

2.3  Thermodynamic  Parameters 

Temperature  Differences.  Figure  8  shows  the  3-5°  radial  belt  aver¬ 
age  of  the  vertical  profile  of  the  temperature  difference  obtained  by 
subtracting  the  non- intensifying  temperature  profile  from  the  intensi¬ 
fying  profile  at  both  the  LF  -42  and  LF  -18  time  periods.  For  both 
time  periods  the  intensifying  set  is  slightly  warmer  through  the  tropo¬ 
sphere  with  the  largest  differences  occurring  in  the  lowest  levels.  At 
the  tropopause  level  the  intensifying  set  is  colder  by  one  to  two 
degrees.  These  differences  could  be  due  to  the  different  average  inten¬ 
sities  between  the  data  sets  but  investigation  of  the  temperature  fields 
at  the  lowest  standard  level  (i.e.,  850  mb)  shows  this  is  not  so 

The  850  mb  temperature  fields  for  both  intensity  classes  at  the 
LF  -42  and  LF  -18  time  periods  are  shown  in  Figs.  9  and  10.  Diagram  c 
of  Figs.  9  and  10  represents  the  difference  between  the  two  data  sets 
found  by  subtracting  the  non- intensifying  data  set  from  the  intensifying 
set.  The  largest  differences  occur  over  the  north  and  northwest  qua¬ 
drants  of  the  storm.  In  this  area  the  intensifying  set  averages  two  to 
three  degrees  warmer  for  both  time  periods.  The  analysis  also  shows  a 
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Fig.  8.  3-5°  radial  belt  average  of  the  vertical  profile  of  the 

temperature  difference  obtained  by  subtracting  the  non¬ 
intensifying  temperature  profile  from  the  intensifying  profile 
at  both  the  LF  -42  and  LF  -18  time  periods. 

temperature  trough  over  the  northwest  sector  for  the  non- intensifying 
set  -  a  feature  not  found  on  the  intensifying  maps. 

Moisture  Differences.  Figure  11  shows  the  3-5°  belt  average  of 
specific  humidity  for  the  LF  -18  time  period  and  Fig.  12  is  the  relative 
humidity  profile  for  the  same  time  and  area.  The  specific  and  relative 
humidities  are  both  slightly  larger  for  the  intensifying  data  set  but 
the  difference  is  so  small  that  it  is  hardly  measurable.  Similar 
comparisons  for  the  LF  -42  time  period  also  indicate  little  or  no  dif¬ 
ference  between  the  two  systems.  This  finding  was  also  observed  by 
Nunez  (1981). 

Height  Differences.  Because  of  the  warmer  troposphere  found  in 
the  intensifying  data  set,  the  hydrostatic  pressure  relationships  dic¬ 
tate  that  the  pressure-height  fields  at  upper  levels  should  be  higher 
in  the  intensifying  set.  The  200  mb  height  analysis  is  shown  in  Figs. 


Fig.  11.  The  3-5  radial  belt  average  of  specific  humidity  for  the  LF 
-18  time  period. 


RELATIVE  HUMIDITY  (%) 


Fig.  12.  The  3-5°  radial  belt  average  of  relative  humidity  for  the  LF 
-18  time  period. 
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13  and  14.  At  the  LF  -42  time  period  (Fig.  13)  the  anticyclone  center 
is  established  in  both  sets  but  the  height  gradient  to  the  north  is 
much  stronger  in  the  non-intensifying  set.  Diagram  c  of  Fig.  13  (ob¬ 
tained  by  subtracting  the  non-intensifying  composite  from  the  intensi¬ 
fying  composite)  shows  a  large  difference  over  the  west,  northwest,  and 
northern  octants.  At  the  LF  -18  time  period  (Fig.  14),  there  is  very 
little  change  in  the  intensifying  set.  The  anticyclone  center  found  on 
the  previous  non- intensifying  figure  has  diminished  and  the  presence  of 
a  trough  impingement  to  the  northwest  is  more  evident.  Again  the  dif¬ 
ferencing  diagram  (Fig.  14c)  shows  large  height  differences  over  the 
west  through  northern  octants. 

Negligible  differences  were  found  in  all  other  thermodynamic  para¬ 
meters. 

2.4  Dynamic  Differences 

Total  Wind.  The  total  winds  are  depicted  in  the  plan  views  shown 
in  Figs.  15-18  for  the  900  and  200  mb  levels  and  for  both  time  periods. 
Differences  in  the  total  winds  at  the  900  mb  level  are  small  for  both 
time  periods.  But  at  200  mb  one  can  see  that  the  westerlies  are  much 
stronger  and  extend  further  south  in  the  non-intensifying  set.  There 
also  is  an  apparent  wind  maximum  between  3-11°  on  the  north  side.  The 
anticyclone  is  not  well  developed  in  the  non- intensifying  set  and  is 
displaced  about  4°  to  the  southeast  at  the  LF  -42  time  period.  By  the 
LF  -18  time  period  the  anticyclone  becomes  even  more  disorganized  and 
a  quick  streamline  analysis  would  put  the  center  at  about  6°  to  the 
southeast.  As  was  previously  found  with  the  200  mb  height  field,  the 
largest  differences  occur  over  the  west  through  north  octants. 


•  't*  .  V-- 
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Zonal  Winds.  The  zonal  wind  at  200  mb  was  next  investigated  and 
the  results  can  be  seen  in  Figs.  19  and  20.  Again  the  major  differences 
occur  to  the  west  and  north  but  the  zonal  wind  also  shows  some  major 
differences  in  the  southwest  octant  providing  more  evidence  of  a  trough 
intrusion.  The  200  mb  wind  maximum  found  in  the  plan  views  of  the  total 
winds  (Figs.  15b  and  17b)  appears  to  be  more  confined  to  the  5-11°  belt 
on  the  north  side.  This  wind  maximum  has  a  strength  of  20-24  ms'*'  and 
is  a  persistent  feature  between  the  two  time  periods. 

Vertical  Wind  Shear.  Another  feature  which  has  been  found  by  other 
investigators  (Gray,  1968,  1975;  McBride,  1979)  and  which  is  also  appar¬ 
ent  in  this  data  set  is  the  presence  of  weak  vertical  wind  shear  near 
the  center  of  the  intensifying  system.  Figures  21  and  22  show  the  plan 
views  of  the  mean  shear  of  the  zonal  wind,  U200  mb  ~  U900  mb  ^or  t*le  tWO 
time  periods.  At  both  time  periods  the  intensifying  set  maintains  the 
zero  shear  line  over  the  top  of  the  storm  while  the  non- intensifying  set 
has  15-20  ms'*  shear  above  storm  center.  McBride  (19  79)  indicated  the 
need  for  strong  westerly  shear  to  the  poleward  side  and  strong  easterly 


shear  to  the  equatorward  side  for  cyclone  development.  This  phenomenon 
is  present  in  the  intensifying  data  set  at  both  the  LF  -42  and  LF  -18 
time  periods.  The  non- intensifying  data  set  does  not  portray  the 


easterly  shear  until  reaching  large  radii  on  the  equatorward  side. 

The  radial  wind,  tangential  wind,  divergence  and  vertical  motion 
parameters  were  closely  examined  for  both  data  sets.  The  differences 
noted  were  that  the  intensifying  systems  maintained  an  inflow  through 
a  deeper  layer  resulting  in  more  convergence  and  stronger  vertical 
motion  than  observed  in  the  non-intensifying  systems.  Because  of  the 
difficulty  in  calculating  the  parameters  on  an  individual  case  basis, 
these  parameters  will  not  be  considered  for  further  investigation. 
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Fig.  15.  Plan  views  of  the  200  mb  total  winds  for  the  intensifying  data 
set  (diagram  a)  and  the  non- intensifying  data  set  (diagram  b) 
for  the  LF  -42  time  period.  Wind  barbs  are  in  knots  where 
a  full  barb  =  10  kts. 


Fig.  16.  Same  as  Fig.  15  but  for  the  900  mb  level. 
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Fig.  17.  Plan  views  of  the  200  mb  total  winds  for  the  intensifying  data 
set  (diagram  a)  and  the  non-intensifying  data  set  (diagram  b) 
for  the  LF  -18  time  period.  Wind  barbs  are  in  knots  where  a 
full  barb  =  10  kts. 
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Fig.  18.  Same  as  Fig.  17  but  for  the  900  mb  level. 


(c) 


Fig.  19.  200  mb  zonal  wind  field  for  the  LF  -42  time  period  for  the 

intensifying  data  set  (diagram  a),  the  non- intensifying  data 
set  (diagram  b)  and  the  difference  between  the  two  (diagram  c) 
found  by  subtracting  the  intensifying  data  set  from  the  non¬ 
intensifying  data  set.  In  diagram  c  positive  numbers  indicate 
that  the  zonal  wind  is  stronger  in  the  non-intensifying  data 
set.  Units  are  m  s”  . 
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2.5  Summary  of  Composite  Results 

The  composite  results  show  that  large  differences  occur  in  specific 
meteorological  parameters  between  intensifying  and  non-intensifying 
tropical  cyclone  systems.  It  appears  that  these  differences  are  a  result 
of  large  scale  baroclinic  interaction  since  most  of  the  differences  were 
found  on  the  poleward  side  of  the  storm.  The  major  results  of  the  com¬ 
posite  study  are  listed  below: 

1)  The  large  scale  environment  of  the  intensifying  cyclone  is 
warmer  through  the  depth  of  the  troposphere  with  the  largest 
differences  found  below  700  mb.  The  horizontal  plan  views  at 
850  mb  indicate  that  most  of  the  temperature  difference  occurs 
over  the  west  through  north  octants. 

2)  Moisture  analysis  is  not  a  good  indicator  of  intensification. 

3)  The  largest  height  differences  occur  at  the  upper  levels  over 
the  west  through  north  octants.  The  horizontal  height  gradients 
are  also  strongest  to  the  north  and  northwest. 

4)  The  westerlies  extend  further  south  and  are  much  stronger  in 
the  non- intensifying  set.  A  small  wind  maximum  occurs  between 
5  and  11  on  the  poleward  side  of  the  non- intensifying  system. 

5)  The  differences  in  the  zonal  wind  are  largest  over  the  south¬ 
west  to  north  octants.  This  yields  evidence  of  a  baroclinic 
trough  impinging  upon  the  non- intensifying  system  from  the 
west. 

6)  An  area  of  zero  vertical  zonal  wind  shear  U20O  -  U^qq 
exists  over  the  center  of  the  intensifying  set  while  15-20 

m  s-l  shear  is  found  over  the  non- intensifying  storm  center. 


3.  DEVELOPMENT  OF  THE  INTENSITY  CHANGE  FORECASTING  SCHEME 


The  results  of  the  composite  study  indicated  important  parameter 
differences  on  the  northwest  side  of  the  cyclone  systems.  Therefore, 
this  study  chose  to  concentrate  on  the  region  between  5  and  11°  over 
octants  8,  1  and  2.  For  cyclones  moving  through  the  Gulf  of  Mexico, 
this  region  gives  the  best  data  coverage  (see  Fig.  23). 

Results  derived  from  the  composite  analysis  of  the  previous  section 
were  applied  to  individual  cases  of  the  composite  data  sets.  This  was 
done  to  see  how  well  the  composite  information  could  be  used  on  an  in¬ 
dividual  case  basis.  Additionally,  it  gives  us  a  quantitative  measure 
for  how  applicable  a  forecasting  scheme  might  be. 

For  simplicity  only  the  standard  weather  charts  analyzed  at  the 
850,  500,  and  200  mb  levels  were  utilized.  These  charts  are  readily 
available  to  forecasters  in  any  weather  station  and  can  be  used  to  make 
the  rapid  computations  necessary  in  an  operational  environment.  Also, 
by  selecting  these  levels  we  can  supplement  the  850  and  200  mb  charts 
with  satellite  derived  winds  in  the  data  sparse  region. 

3.1  Selection  of  Parameters 

Temperature.  The  composite  results  indicated  that  the  850  mb  tem¬ 
perature  field  was  noticeably  different  over  the  northwest  octant  of 
the  intensifying  vs.  non-intensifying  storm  system.  The  850  mb  tempera¬ 
ture  was  selected  as  one  of  the  parameters  and  an  average  value  was 
computed  over  the  octants  8,  1  and  2  within  the  5-11°  radial  belt 
through  linear  averaging  by  2°  radial  belts.  Applying  this  procedure 
to  the  two  data  sets  at  both  time  periods  yielded  an  average  850  mb 
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southeast  to  west  Gulf  regions. 

temperature  of  17.5°C  for  the  intensifying  data  set  vs.  15.5°C  for  the 
non- intensifying  data  set  or  a  difference  of  2°C. 

Height.  The  200  mb  height  differences  between  the  intensifying 
and  non- intensifying  data  sets  computed  over  the  5-11°  radial  belt  over 
octants  8,  1  and  2  was  about  60  m  for  both  the  LF  -42  and  LF  -18  time 
periods.  Heights  were  12410  m  for  the  intensifying  set  and  12350  m  for 
the  non-intensifying  set.  The  5-11°  radius  horizontal  height  gradients 
were  analyzed  but  the  differences  between  the  data  sets  was  only  about 
12  m  per  6°  radius. 

Zonal  Wind.  Large  differences  were  noted  earlier  in  the  zonal 
wind  fields.  Vertical  profiles  of  the  zonal  wind  were  first  made  for 
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the  5-7°  radial  belt  over  the  northwest  sector,  octants  1,  2  and  3  (Fig 
24).  Very  large  differences  occur,  particularly  in  the  strength  of  the 
westerlies  in  the  upper  levels  and  in  the  vertical  shear.  But  for 
storms  in  the  Gulf  of  Mexico  we  encountered  the  problem  of  data  defi¬ 
ciencies  over  octant  3.  Therefore,  a  similar  profile  was  made  for  the 
5-7^  radial  belt  over  octants  8,  1  and  2.  Results  are  shown  in  Fig.  25 
Although  the  extrema  are  not  as  large,  they  are  still  very  striking, 
especially  at  200  mb  where  the  differences  are  more  than  10  m  s  ^  be¬ 
tween  systems. 

Vertical  Wind  Shear.  Finally,  to  be  consistent  with  the  methods 
for  computing  the  temperature  and  height  values  and  also  to  smooth  out 
possible  eddies  in  the  6°  radial  belt,  a  vertical  profile  of  the  zonal 
wind  was  made  for  the  5-11°  belt  over  octants  8,  1  and  2  as  shown  in 
Fig.  26.  The  same  differences  in  the  zonal  wind  still  occur  although 
they  are  not  as  large.  By  expanding  the  areal  coverage  of  the  measure- 
roent  we  can  minimize  the  data  noise  and  get  a  more  stable  representa¬ 
tion  of  the  parameter  differences.  Therefore,  we  will  average  all 
parameters  over  the  5-11°  radial  belt.  This  improves  the  data  quality 
without  greatly  reducing  the  differences  between  the  two  systems. 

The  vertical  wind  shears  computed  over  the  5-11°  radial  belt  over 
octants  8,  1,  and  2  between  850-500  mb  and  850-200  mb  are  also  very 
different  (Fig.  26).  The  850-500  mb  vertical  wind  shear  is  only  about 
3  m  s  ^  for  the  intensifying  set  vs.  6  m  s  *  for  the  non- intensifying 
set.  The  850-200  mb  vertical  shear  is  about  12  m  s  ^  for  the  inten¬ 
sifying  set  vs.  21-25  m  s  ^  for  the  non-intensifying  data  set.  The 
differences  are  even  larger  if  we  compare  Figs.  24  and  25. 
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24.  The  vertical  profile  of  the  zonal  wind  for  the  6°  radial  belt 
averaged  over  octants  1,  2  and  3  for  the  LF  -42  time  period 
(left  diagram)  and  the  LF  -18  time  period  (right  diagram). 


Fig.  25.  The  vertical  profile  of  the  zonal  wind  for  the  6°  radial  belt 
averaged  over  octants  8,  1  and  2  for  the  LF  -42  time  period 
(left  diagram)  and  the  LF  -18  time  period  (right  diagram). 


Fig.  26.  The  vertical  profile  of  the  zonal  wind  averaged  over  the 
5-11  radial  belt  for  octants  8,  1  and  2  for  the  LF  -42 
time  (left  diagram)  and  the  LF  -18  time  period  (right  diagram). 


Summary.  Further  analysis  of  the  composite  data  sets  on  the  north¬ 
ern  side  continues  to  produce  large  differences  between  the  two  data 
sets.  But  the  unanswered  question  is  how  much  variability  exists  in 
the  parameters  of  the  individual  cyclone  cases  which  make  up  the  compo¬ 
site  study.  Furthermore,  how  well  can  we  use  these  results  to  produce 
skillful  forecasts?  To  answer  these  questions,  the  parameters  discussed 
in  this  section  will  be  measured  for  each  individual  cyclone  system  and 
compared  against  the  mean  to  see  if  we  can,  in  fact,  produce  a  reliable 
forecasting  scheme. 

3.2  Application  of  Parameters  to  Individual  Cases 

Special  computer  programming  techniques  developed  by  Edwin  Buzzell 
of  the  CSU  tropical  storm  group  were  employed  to  retrieve  the  data  for 
each  of  the  individual  storms  which  made  up  the  composites.  Objective 
procedures  were  used  to  calculate  all  parameters.  In  data-sparse  regions 
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smoothing  procedures  were  used  to  determine  the  parameter  values.  A 
complete  listing  of  the  computations  can  be  found  in  Tables  7-10  of 
Appendix  A.  A  summary  of  those  tabular  values  can  be  seen  in  Figs.  27 
and  28  and  will  be  discussed  in  the  following  paragraphs. 

Temperat ure.  The  average  5-11°  radius  850  mb  temperature  for  all 
cyclones  ranged  from  13. 5°C  to  21°C.  The  value  of  15.5°C  was  chosen  as 
a  threshold  value  upon  which  to  base  intensification.  Only  9%  of  the 
intensifying  systems  (2  of  23  cyclones  averaged  over  the  LF  -42  and  LF 
-18  time  periods)  had  an  850  mb  temperature  colder  than  15.5°C.  Of  the 
non- intensifying  systems,  50%  (6  of  12  storms)  have  850  mb  temperatures 
warmer  than  15.5°C  but  had  other  parameters  which  did  not  favor  inten¬ 
sification.  Of  all  cyclones  averaged  over  the  two  time  periods,  77% 
could  have  been  identified  as  intensifying  or  non- intensifying  simply 
from  their  respective  5-11°  radius  850  mb  temperature  fields  being 
warmer  or  colder  than  15.5°C. 

Height.  The  200  mb  height  field  averaged  over  the  5-11°  radial 
belt  over  octants  8,  1  and  2  also  displays  a  good  delineation  between 
systems.  The  value  of  12380  m  provided  a  base  value  upon  which  74%  and 
80%  of  the  storms  at  the  LF  -42  and  LF  -18  time  periods  could  have  been 
characterized  as  intensifying  or  not.  Most  of  the  remaining  26-20% 
was  made  up  of  intensifying  systems  with  height  values  <  12380  m.  In 
no  case  did  a  storm  intensify  with  a  height  field  <  12320  m  and  con¬ 
versely  only  one  storm  did  not  intensify  when  the  height  field  was 
>  12420  m.  Furthermore,  at  the  LF  -18  time  period  as  the  cyclone 
approached  landfall  and  further  baroclinic  interaction,  only  one  cyclone 
did  not  intensify  with  a  height  field  >  12380  m. 


Fig.  27.  Tabular  values  of  each  of  the  five  intensity  change  parameters  as  they  are  distributed  in  their 
respective  scatter  graph  for  the  LF  -42  time  period.  The  0  line  represents  the  threshold  value 
of  each  parameter  as  determined  from  the  scatter  of  the  data.  The  abscissa  on  the  bottom  of  the 
graph  is  used  in  the  computation  of  the  intensity  change  parameter. 


as  Fig.  27  but  for  the  LF  -18  time  period. 
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Zonal  Wind.  The  individual  case  measurements  of  the  impinging 
zonal  westerlies  upon  the  storm  system  also  verifies  the  results  of  the 
composites.  The  negative  of  the  zonal  wind  (-  I^qq)  is  plotted  in  Figs. 
27  and  28  and  the  established  threshold  is  -10  m  s  \  The  results  are 
best  at  the  LF  -18  time  period  where  80%  of  the  cyclones  are  properly 
classified  as  intensifying  or  not.  The  majority  of  intensifying 
cyclones  which  were  incorrectly  classified  had  zonal  winds  stronger  than 
the  threshold  value.  These  storms  were  primarily  large  strong  storms 
with  a  well  established  anticyclone.  The  effects  of  the  strong  anti¬ 
cyclone  at  200  mb  was  to  add  to  the  strength  of  the  westerlies  in  the 
6°  belt  on  the  northern  octants. 

Vertical  Wind  Shear.  Similar  results  were  also  found  in  the 
vertical  zonal  wind  shear  calculations.  A  threshold  value  of  -16  m  s  ^ 
for  the  850  mb  minus  200  mb  zonal  wind  was  established  from  the  scatter 
of  the  data  as  depicted  in  Figs.  27  and  28.  Here  again  the  intensifying 
cyclones  with  large  vertical  wind  shear  were  primarily  large  cyclones 
with  a  strong  upper  level  anticyclone.  However,  most  of  these  cyclones 
had  a  compensating  weak  vertical  wind  shear  from  850  to  500  mb. 

To  summarize  Figs.  27  and  28  we  find  that  a  large  variability  exists 
in  each  of  the  parameters  but  a  distinct  division  is  evident  between  the 
two  data  sets.  An  attempt  to  maximize  the  positive  benefits  of  each 
parameter  was  made  by  adjusting  the  scale  of  each  parameter  based  on 
its  importance  as  determined  from  statistical  methods  shown  in  section  5. 

3.3  Designation  of  the  Forecasting  Parameter 

An  intensification  parameter  (IP)  can  be  determined  from  a  normal¬ 
ized  sum  of  each  parameter  as  determined  from  the  scale  on  the  bottom 
of  Figs.  27  and  28  or 
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IP  =  H200  +  T850  +  U200)  +  (U850  ~  U200)  +  (U850  ~  U50(P  (2) 

This  IP  parameter  was  designed  to  imply  the  nature  of  a  storm's  poten¬ 
tial  intensification  but  not  necessarily  the  degree  of  that  intensifica¬ 
tion.  For  example,  if  IP  <  0  we  would  expect  this  storm  to  be  a 
non- intensifying  system;  if  IP  >  0  the  storm  should  intensify  prior  to 
landfall.  There  is  a  climatological  variability  in  IP  through  the 
storm  season  and  an  adjustment  for  this  variability  is  made  to  compen¬ 
sate  for  the  imbalance.  Discussion  of  the  climatic  variability  is  given 
in  section  4.2  and  details  of  the  correction  are  discussed  in  Appendix 
B. 

IP  applied  at  the  LF  -42  time  period  correctly  forecasts  the  inten¬ 
sity  change  for  71%  of  the  35  storms.  At  the  LF  -18  time  period  the 
skill  improves  to  77%.  IP  was  plotted  against  the  next  24  hour  wind 
speed  change  for  each  storm  and  the  results  can  be  seen  in  Figs.  29  and 
30.  The  statistical  implications  of  these  figures  will  be  treated  in 
section  5  but  the  two  data  sets  appear  significantly  different.  The 
exceptions  of  particular  note  in  the  figures  are  the  open  circles 
located  to  the  left  of  the  zero  line.  These  data  points  represent  the 
following  storms:  Cindy  1963,  Isbell  1964,  Abby  1968,  Edith  1971, 

Carmen  1974,  and  Eloise  1975.  A  careful  analysis  of  the  weather  charts 
for  these  storms  revealed  that  they  all  moved  into  and  intensified  in  a 
baroclinic  region.  This  is  precisely  the  kind  of  interaction  which 
typically  prevented  the  non-intensifying  systems  from  intensifying. 
Perhaps  these  storms  had  different  inner  core  characteristics  and  the 
baroclinic  interaction  in  some  way  complemented  the  dynamics  of  the 
storm.  It's  possible  that  a  measurement  of  the  imposing  baroclinic 
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Fig.  29.  Intensity  change  parameter  (IP)  vs.  the  next  24  hour 

sustained  wind  speed  change  for  the  LF  -42  time  period. 
Circled  data  points  are  the  central  values  of  the  respective 
set. 
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field  alone  is  not  enough  and  that  there  is  a  combination  of  several 
forces  which  governs  intensity  change.  Careful  budget  studies  (parti¬ 
cularly  angular  momentum)  may  yield  a  better  answer  to  each  storms’ 
intensification  trends. 

3.4  Rate  of  Change  of  IP 

Another  important  feature  of  the  intensity  change  parameter  is 
its  rate  of  change.  Preliminary  results  of  case  studies  show  that  when 
the  intensity  parameter  registers  a  large  increase  with  time,  the  storm 
typically  responds  with  a  strong  intensification  within  6  to  12  hours. 
Conversely,  a  sharp  decrease  in  IP  is  usually  followed  by  a  steady  state 
or  filling  storm.  Therefore,  the  initial  value  of  IP  should  not  be 
the  sole  determinator  in  forecasting  storm  intensity  but  the  trend  of 
the  IP  parameter  should  also  be  taken  into  consideration. 
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4.  APPLICATION  OF  IP  TO  INDEPENDENT  DATA 

It  would  be  unfair  to  make  general  conclusions  about  the  reliability 
of  the  forecasting  parameter  based  on  only  35  cases.  Therefore,  the 
author  sought  independent  data  on  which  to  further  test  the  parameter. 
Seven  cyclones  from  the  storm  seasons  of  1978  and  1979  were  examined. 
Secondly,  all  tropical  depressions  for  the  period  1967-1979  (23  cases) 
were  examined.  Then,  the  climatological  and  diurnal  variabilities  on 
the  parameters  were  examined.  Included  is  a  discussion  of  the  adjust¬ 
ments  made  to  IP  to  allow  for  the  variabilities.  Finally,  the  scheme 
was  applied  to  23  independent  intensity  change  composite  data  sets  of 
the  west  Atlantic  and  west  Pacific  available  at  CSU.  A  listing  of  all 
additional  cyclones  and  the  analysis  of  their  various  parameters  is 
given  in  Tables  11-14  of  Appendix  A. 

4.1  Results  of  Independent  Storm  Analysis 

All  independent  cyclone  track  positions  were  obtained  from  the 
annual  best  track  summaries.  Tropical  depressions  were  first  tracked 
in  1967  but  it  wasn't  until  1974  that  they  were  assigned  numbers.  In¬ 
tensity  changes  for  depressions  are  not  included  in  the  summaries  until 
they  reach  a  minimum  tropical  storm  strength.  Hence,  the  author  had  to 
make  a  general  assumption  that  all  depressions  were  non- intensifying 
systems.  This  is  not  necessarily  a  good  assumption  since  depressions 
have  been  known  to  form  in  weak  wave-trough  systems,  intensify  to 
15  is  s  ^  within  24  hours  and  then  make  landfall  a  short  time  later  — 
all  criteria  which  would  classify  the  system  as  intensifying.  But  for 
the  purposes  of  this  study  all  depressions  will  be  treated  equally  with 
the  exception  of  Depression  #9  of  1975.  This  depression,  as  noted  by 
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N.  Frank  (1975),  was  one  which  formed  quickly  and  caused  millions  of 
dollars  of  damage  in  flooding  as  it  moved  inland.  Although  exact  wind 
speed  change  was  not  available  for  this  depression,  indications  were 
that  it  was  an  intensifying  system. 

The  parameters  were  analyzed  for  all  independent  cyclones  using  the 
individual  weather  charts  following  the  computational  procedures  for 
intensity  change  as  discussed  in  Appendix  B.  The  results  of  this 
analysis  are  presented  in  Figs.  31  and  32.  Since  data  on  the  next  24 
hour  wind  speed  change  was  unknown  for  the  depressions,  it  was  assumed 
that  these  wind  speed  changes  were  zero.  Even  with  the  assumption  that 
all  depressions  are  non- intensifying  systems  we  could  have  forecast  the 
correct  intensification  trend  for  15  of  20  cyclones  (75%)  at  the  LF  -42 
time  period  and  24  of  31  cyclones  (77%)  at  the  LF  -18  time  period. 

These  results  are  comparable  to  those  obtained  from  the  original  data 
sets.  The  big  difference  in  examining  the  independent  cyclone  data  is 
that  the  false  alarm  rate  has  increased  (more  cyclones  forecast  to  in¬ 
tensify  than  actually  did)  but  the  rate  of  missing  an  intensity  increase 
has  been  lowered. 

Tables  7-14  of  Appendix  A  list  the  data  for  each  cyclone  case.  To 
give  a  better  understanding  of  the  variability  of  the  IP  parameter  for 
each  data  set  as  listed  in  the  tables,  Fig.  33  was  constructed  to  show 
how  the  data  is  distributed  about  the  median  value.  The  solid  line 
represents  the  range  over  which  approximately  50%  of  all  cases  fall. 

The  dashed  line  gives  the  range  in  which  about  70%  of  all  cases  fall. 
Because  of  the  small  sample  size  (5  cases)  of  the  independent  inten¬ 
sifying  data  set  at  the  LF  -42  time  period,  the  solid  line  for  this  data 
set  represents  a  60%  range.  Figure  33  shows  no  overlap  between 
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31.  Intensity  change  parameter  (IP)  versus  the  next  24  hour 
sustained  wind  speed  change  for  all  independent  storms  at 
the  LF  -42  time  period.  Circled  data  points  are  the  central 
values  of  the  respective  data  sets. 
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Fig.  32.  Same  as  Fig.  31  but  for  the  LF  -18  time  period. 
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DATA  SET 

SAMPLE  : 

1. 

Intensifying  LF  -42 

23 

2. 

Non- intensifying  LF  -42 

12 

3. 

Intensifying  LF  -18 

23 

4. 

Non- intensifying  LF  -18 

12 

5. 

Independent  intensifying  LF  -42 

5 

6. 

Independent  non- intensifying  LF  -42 

15 

7. 

Independent  intensifying  LF  -18 

7 

8. 

Independent  non- intensifying  LF  -18 

24 

9. 

Average  of  all  intensifying  cases 

58 

10. 

Average  of  all  non-intensifying  cases 

63 

Fig.  33.  Variability  of  the  IP  parameter  for  individual  cyclone  cases 
as  distributed  about  the  median  value  of  each  data  set.  The 
solid  line  represents  the  range  over  which  approximately  50% 
of  all  cases  fall.  The  dashed  line  gives  the  range  in  which 
about  70%  of  all  cases  fall.  An  exception  is  data  set  5  which 
represents  a  60%  range  because  of  the  small  sample  size.  The 
number  beside  the  open/closed  circles  represents  the  data 
sets  described  above. 


intensifying  and  non-intensifying  cases  at  the  50%  level.  However,  at 
the  70%  level  some  overlap  occurs  in  nearly  all  data  sets.  This  is  the 


gray  region  in  which  the  forecast  errors  of  this  scheme  originate. 


4.2  Climatological  and  Diurnal  Variation  of  the  Parameters 

The  monthly  distribution  of  intensifying  vs.  non-intensifying  cases 
(Table  2  -  page  15)  gives  indications  that  the  differences  found  in  the 
composite  study  may  be  due  to  climatological  differences  found  by  com¬ 
paring  early  and  late  storms  to  those  which  occur  during  the  peak  of  the 
storm  season.  However,  after  the  addition  of  the  independent  cyclone 
data  almost  no  such  seasonal  variability  exists.  Table  3  shows  nearly 
an  equal  distribution  of  cyclones  by  month  for  both  the  intensifying 
and  non- intensifying  cyclone  systems. 

Seasonal  climatological  variations  do  occur  in  all  parameters. 
Therefore,  the  IP  parameter  must  also  have  a  seasonal  variation  which 
must  be  compensated  for  in  the  forecast  scheme.  A  ten  year  (1958-1967) 
monthly  mean  over  18  rawinsonde  stations  located  in  the  Gulf  coastal 
area  was  computed  for  H^qq  and  Tg,.Q.  The  seasonal  variation  of  these 
parameters  is  shown  in  Figs.  34  and  35.  The  dashed  line  in  the  figure 
is  the  threshold  value  established  for  the  forecast  scheme.  Monthly 
mean  values  of  the  wind  parameters  (Ug^Q»  U5qo  and  U200^  were  evaluated 
by  averaging  5  Gulf  coastal  rawinsonde  stations  over  a  5  year  period 
(1960-1964).  From  these  climatological  mean  values  of  all  parameters, 
a  monthly  mean  value  of  the  IP  parameter  was  calculated.  The  variations 
of  the  IP  parameter  is  shown  by  the  heavy  solid  line  in  Fig.  36.  All 
cyclone  cases  from  this  study  were  subdivided  into  15  day  (half-month) 
periods  and  average  IP  values  were  computed  for  the  intensifying  cases, 
the  non- intensifying  cases  and  the  average  of  all  cases  as  shown  in 
Fig.  35.  The  seasonal  variations  of  all  individual  cases  closely 
parallels  the  curve  determined  from  climatology.  A  compensation  curve 
to  allow  for  these  seasonal  variations  in  the  IP  parameter  is  shown  in 
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TABLE  3 

Monthly  occurrence  of  all  cyclones  which  make  up  the  study.  The  table 
includes  both  the  individual  cyclones  from  the  composite  study  and  the 
cyclones  from  the  independent  sample. 


MAY  JUN  JUL  AUG  SEP  OCT  TOTAL 


INTENSIFYING 

0 

□ 

D 

6 

14 

2 

30 

NON-INTENSIFYING 

2 

5 

5 

6 

12 

6 

36 

Fig.  37.  Basically,  the  calculated  value  of  the  IP  is  increased  during 
the  early  and  late  portion  of  the  storm  season  and  decreased  during  the 
peak  of  the  season.  Exact  details  of  the  corrections  are  given  in 
Appendix  B. 

Diurnal  variations  in  Tg^  were  noted  when  making  the  computations 
for  the  independent  cyclone  systems.  The  variations  were  particularly 
significant  for  cyclones  centered  in  the  western  Gulf  of  Mexico.  For 
cyclones  in  this  region,  the  5-11°  belt  of  octant  2  is  located  over 
central  and  west  Texas.  This  part  of  Texas  is  slightly  elevated  and 
experiences  strong  surface  heating.  The  surface  heating  is  reflected 
at  850  mb  where  temperatures  as  high  as  27°C  were  observed  at  00Z  in 
contrast  to  an  average  temperature  of  17°C  at  12Z.  For  storms  in  this 
area,  only  the  12Z  temperature  should  be  used.  Large  variations  in 
temperature  were  not  observed  in  other  areas. 

Diurnal  variations  in  could  not  be  detected. 

4.3  Comparative  Composite  Analysis 

Accepting  the  hypothesis  that  poleward  baroclinic  interaction 
causes  storm  weakening,  it  is  to  be  expected  that  trends  similar  to 
those  shown  for  the  Gulf  of  Mexico  should  also  occur  in  other  ocean 
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Fig.  34.  Ten  year  (1958-1967)  monthly  mean  of  the  200  mb  height 

averaged  over  18  rawinsonde  stations  through  the  Gulf  coastal 
area.  The  dashed  line  is  the  threshold  value  established  for 
this  study. 


Fig.  35.  Same  as  Fig.  34  but  for  the  850  mb  temperature. 


5 


Fig.  36.  Mean  monthly  values  of  IP  as  computed  from  climatology  (heavy 

solid  curve) ,  all  individual  intensifying  cases  (dotted  curve) , 
all  individual  non-intensifying  cases  (dashed  curve)  and  an 
average  of  both  intensifying  and  non-intensifying  cases 
(dash-dotted  curve) 


Fig.  37.  Seasonal  correction  curve  to  be  applied  in  the  computation  of 
IP. 
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basins.  The  author  examined  the  comparative  composite  studies  made  by 
previous  investigators  at  CSU  —  Frank  and  Gray  (1976),  Zehr  (1976), 
McBride  (1979)  and  Ntinez  (1981).  Only  those  studies  which  contained  a 
number  of  storms  poleward  of  15°  latitude  were  considered.  Climatolo- 
gically,  baroclinic  troughs  rarely  penetrate  equatorward  of  15° 
latitude  during  the  storm  season.  See  Appendix  C  for  a  detailed  descrip¬ 
tion  of  each  composite  data  set. 

Analysis  of  these  comparative  composite  studies  (Table  4)  yields 
results  similar  to  those  found  in  the  Gulf  of  Mexico.  In  all  filling 
composite  data  sets  the  intensity  change  parameter  is  significantly 
smaller  than  for  the  corresponding  deepening  set.  This  difference  can 
easily  be  attributed  to  baroclinic  interaction  by  examining  the  indivi¬ 
dual  parameters  of  Table  4.  For  example,  for  the  non-intensifying  cases 
the  westerlies  are  much  stronger  and  are  closer  to  the  storm  center; 
the  vertical  wind  shears  are  much  stronger;  the  low  level  temperature 
fields  are  colder;  and  the  upper  level  height  fields  are  lower. 

Minor  modifications  had  to  be  made  to  the  procedures  for  computing 
IP  in  the  Pacific  basin.  Climatologically ,  the  850  mb  temperature  field 
is  warmer  and  the  200  mb  height  field  is  higher  in  the  Pacific  than  the 
Atlantic.  Therefore,  new  threshold  values  had  to  be  established  for 
these  parameters.  IP  was  computed  for  all  data  sets  and  the  results 
can  be  seen  in  Fig.  38.  It  is  most  important  to  compare  the  distance 
separating  each  of  the  data  sets  rather  than  the  relative  placement  on 
the  graph  as  in  the  AIP  column  of  Table  4.  Further  study  in  the  Pacific 
basin  needs  to  be  done  on  individual  cases  to  better  refine  the  thres¬ 
hold  values  but  the  preliminary  results  of  Fig.  38  show  great  promise 
for  implementation  of  a  quantitative  system  for  measuring  intensity 
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Fig. 

38.  Graph  of  IP  for  comparative  intensity  change  rawinsonde  com- 

posite  data  sets.  Numbers  denote  each  composite  data  set 


comparison.  The  line  between  the  data  points  shows  the  dif¬ 
ference  between  the  intensifying  and  non-intensifying  data 
sets.  See  Appendix  C  for  a  detailed  description  of  all  data 
sets. 
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change.  The  climatological  and  diurnal  effects  must  also  be  examined 
more  carefully  in  this  basin  before  a  scheme  such  as  this  can  be  used. 

The  forecasting  scheme  was  also  applied  to  a  set  of  rapidly  deep¬ 
ening  typhoons  (deepening  rate  >  42  mb  d  from  the  northwest  Pacific 
which  was  studied  by  Holliday  and  Thompson  (1979).  IP  was  determined 
from  rawinsonde  composites  prepared  for  these  storms  for  three  time 
periods  (RD-24,  RD-0,  RD+24)  as  described  in  Appendix  C.  Computations 
of  IP  made  for  these  Holliday  and  Thompson  cases  indicated  that  inten¬ 
sification  should  occur.  Close  analysis  of  the  parameters  which  deter¬ 
mine  IP  (Table  4)  revealed  that  a  large  part  of  the  increase  in  IP 
between  time  periods  for  this  data  set  was  due  to  changes  in  the  200  mb 
height  field.  It  appeared  that  the  storms  were  moving  into  an  area  of 
higher  200  mb  heights  (>  12470  m)  as  the  rapid  deepening  was  occurring. 
Indications  are  that  it  may  be  just  as  important  to  watch  for  changes 
in  the  individual  parameters  as  well  as  the  total  change  in  IP. 

Discussion.  This . section  has  provided  more  verification  that 
baroclinic  interaction  does  in  fact  modify  a  storm's  intensity  and  more 
importantly  that  the  interaction  can  be  measured  in  a  quantitative  sense. 
The  skill  of  making  forecasts  using  this  quantitative  scheme  has  been 
shown  to  be  >  70%  correct.  The  author  is  hopeful  that  through  further 


refinements  this  skill  can  be  increased. 


5.  STATISTICAL  ANALYSIS 


The  data  as  presented  in  the  study  appear  to  be  highly  significant. 
To  determine  the  level  of  significance  of  these  results,  all  data  were 
examined  collectively  for  the  time  periods  LF  -42  and  LF  -18.  The  sta¬ 
tistical  methods  used  for  these  non-parametric  data  sets  were  contingency 
tables,  prefigurance  -  post  agreement,  skill  scores  and  the  recently 
developed  technique  termed  Multi-Response  Permutation  procedures  (MRPP) . 

A  complete  description  of  the  MRPP  technique  can  be  found  in  Mielke 
et  al.  (1976,  1981)  and  Mielke  (1979). 


5.1  Contingency  Tables 

Contingency  tables  represent  a  means  of  testing  the  discrepancy 
between  observed  results  and  those  expected  under  the  proposed  forecast¬ 
ing  parameter.  Two  by  two  contingency  tables  were  used  to  plot  future 
forecast  intensification  (I)  and  non- intensification  (NI)  against  that 
actually  observed  (Fig.  39).  For  this  statistical  test  any  cyclone 
which  experienced  a  wind  speed  change  >  10  kts  over  the  next  24  hour 
period  was  considered  an  intensifying  system.  Figure  40  shows  the  no¬ 
skill  table  that  would  have  resulted  if  the  forecaster  really  had  no 
skill  and  his  correct  forecasts  were  due  to  random  probability.  This 
no-skill  table  is  formed  from  Fig.  39  by  multiplying  row  totals  by 
column  totals  and  then  dividing  by  the  grand  total.  Using  these  two 
tables  we  can  then  compute  the  chi-square  statistic  by 
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Fig.  39.  Contingency  tables  made  from  all  storm  cases  for  the  LF  -  42 
time  period  (left  diagram)  and  the  LF  -18  time  period  (right 
diagram) . 
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Fig.  40.  No-skill  table  made  from  the  contingency  tables  in  Fig.  39 

by  multiplying  the  row  total  by  the  column  total  and  dividing 
by  the  grand  total.  Left  diagram  is  for  LF  -42  time  period 
and  right  diagram  is  for  the  LF  -18  time  period. 
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where  0^  is  the  observed  counts  (Fig.  39)  and  1L  is  the  no-skill  counts 

2 

(Fig.  40).  We  find  that  for  the  time  period  LF  -42,  x  =12.4  and  at 
2 

LF  -18,  x  =  19.7.  Using  chi-square  tables  with  one  degree  of  freedom, 
we  find  that  for  both  time  periods  the  probability  value  for  chi-square 
is  less  than  0.001  that  a  no-skill  forecast  would  have  led  to  the  given 
verification  tables  (Fig.  39). 

The  usefulness  of  the  predictor  can  also  be  determined  by  forming 
contingency  ratios  (Fig.  41).  These  are  simply  the  ratios  of  the 
numbers  in  the  contingency  tables  (Fig.  39)  to  the  numbers  in  the  cor¬ 
responding  no-skill  tables  (Fig.  40).  Whenever  these  ratios  are  greater 
than  unity,  the  relationship  between  the  predictor  and  outcome  is  much 
better  than  a  chance  one.  As  we  can  see  from  Fig.  41  the  forecasting 
parameter  is  doing  better  than  chance  in  forecasting  both  intensifica¬ 
tion  and  non-intensification. 


5.2  Prefigurance  -  Post  Agreement 


Prefigurance  measures  the  extent 
notice  of  the  occurrence  of  a  certain 


to  which  forecasts  give  advance 
event.  Post  agreement  gives  the 
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Fig.  41.  Contingency  ratios  formed  by  dividing  the  values  in  the  con¬ 
tingency  tables  (Fig.  39)  by  the  values  in  the  no-skill  tables 
(Fig.  40).  Left  diagram  is  for  LF  -42  time  period  and  right 
diagram  is  for  LF  -18  time  period. 
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extent  to  which  subsequent  observations  confirm  the  predictions  when 
a  certain  event  is  forecast.  This  information  can  be  found  by  using 
the  previously  constructed  contingency  tables  (Fig.  39)  to  develop  the 
tables  found  in  Figs.  42  and  43.  The  data  from  these  figures  is  con¬ 
densed  below: 

TABLE  5 

Prefigurance  Post  Agreement 

(Percentage  Accuracy)  (Percentage  Accuracy) 


1 

NI 

I 

NI 

LF 

-42 

80 

63 

65 

79 

LF 

-18 

81 

74 

74 

81 

As  an  example  of  the  use  of  Table  5,  we  might  consider  the  LF  -42  time 
period.  We  see  that  forecasts  of  intensification  were  followed  by  in¬ 
tensification  65%  of  the  time  while  occurrences  of  intensification  were 
indicated  in  advance  80%  of  the  time.  The  trend  in  the  scheme  is  to 
over  forecast  intensification  while  under  forecasting  non- intensif cation. 

5.3  Skill  Scores 

The  information  contained  in  the  contingency  tables  can  be  combined 
into  a  single  index  (S) ,  called  a  skill  score.  It  is  defined  by 


S 


C  -  E 
T  -  E 


where  C  is  the  number  of  correct  forecasts,  T  is  the  total  number  of 
forecasts  and  E  is  the  number  expected  to  be  correct  by  chance.  Skill 
scores  were  computed  for  the  two  time  periods  and  found  to  be  .43  and 
.55  for  the  LF  -42  and  LF  -18  time  periods  respectively.  These  numbers 
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Fig.  42.  Percent  of  time  each  observed  category  was  correctly  forecast 
(Pref igurance) . 
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Fig.  43.  Percent  of  time  each  forecast  event  occurred  for  a  particular 
category  (Post  Agreement) . 
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represent  an  average  correct  forecast  rate  of  75%  which  beats  chance  by 
25%. 


5.4  MRPP 

The  MRPP  technique  which  has  been  developed  at  CSU  (cf.  Mielke  et. 
al. ,  1976,  1981;  Mielke,  1979)  is  a  nonparametric  statistical  methodology 
which  does  not  involve  normality  or  other  assumptions  usually  made  in 
statistical  analyses.  It  allows  the  simultaneous  comparison  of  one  or 
more  variates  of  one  data  set  with  the  same  variates  of  other  data  sets. 
The  procedure  essentially  measures  the  euclidean  distance  between  all 
data  points  giving  the  likelihood  of  one  set  of  data  being  statistically 
the  same  as  an  opposing  data  set.  The  MRPP  technique  was  applied  to  the 
data  of  this  study  in  three  ways: 

1)  All  parameters  were  examined  univariately  comparing  the  inten¬ 
sifying  versus  non- intensifying  set. 

2)  The  three  most  significant  parameters,  H200>  T850  an<^  u200 
(6  belt)  were  combined  in  a  three-way  combination. 

3)  All  five  parameters  were  examined  in  a  five-way  combination. 

The  results  of  the  analysis  are  shown  in  Table  6.  One  can  see  that  the 
vertical  wind  shears  are  the  weakest  parameters  but  when  combined  with 
the  other  parameters  to  form  IP,  the  results  become  extremely  signifi¬ 
cant. 

The  MRPP  statistic  can  also  be  used  to  improve  the  selection  of 
parameters  in  devising  a  forecasting  scheme.  For  example,  in  the  NW 
Pacific,  data  is  less  abundant  particularly  at  the  500  mb  level.  There¬ 
fore,  use  of  the  scheme  without  the  contribution  of  the  850  to  500  mb 
shear  should  also  give  good  results.  In  fact  the  MRPP  results  indicate 
that  the  200  mb  height  field  alone  may  be  sufficient  to  produce  a 


reasonable  forecast. 


5  Way  combination  consists  of  all  five  parameters. 


3  Way  combination  consists  of  H Toc„  and  (6  belt) 


6.  CONCLUSION 


This  study  was  undertaken  to  compare  intensifying  and 
non-intensifying  tropical  storm  systems.  The  original  goals  were  to 
detect  large  scale  meteorological  differences  between  the  two  types  of 
storms  merely  by  examining  the  weather  charts.  But  hours  of  tedious 
effort  spent  analyzing  the  charts  could  only  give  subjective  indications 
of  possible  baroclinic  interaction  with  the  cyclone  systems.  Therefore, 
the  rawinsonde  compositing  approach  was  used  to  look  for  systematic  dif¬ 
ferences  in  an  objective  manner.  The  results  obtained  were  significant. 
Of  particular  interest  were  the  differences  found  in  the  850  mb  tempera¬ 
ture  field,  the  200  mb  height  field,  the  200  mb  zonal  wind  and  vertical 
wind  shears  of  the  zonal  wind. 

Based  on  the  results  obtained  from  the  composite  study,  the  indi¬ 
vidual  cyclones  were  again  analyzed  for  the  specific  parameters.  Some 
variability  was  inherent  in  the  results,  but  a  distinct  demarcation 
between  the  two  data  sets  was  evident.  Special  combinations  of  the 
parameters  led  to  a  quantitative  scheme  for  forecasting  tropical  cyclone 
intensity  change. 

The  forecasting  scheme  held  up  well  when  tested  against  independ¬ 
ent  data.  The  scheme  was  able  to  forecast  intensity  change  71-77% 
correct  for  all  time  periods  but  was  at  its  best  as  the  cyclones  were 
within  24  hours  of  landfall.  It  was  found  that  trends  in  the  para¬ 
meters  should  also  be  carefully  watched  as  indicators  of  change. 

Through  the  use  of  previous  comparative  rawinsonde  composite  studies, 
we  were  able  to  show  that  the  results  of  this  study  appear  to  be 
applicable  in  other  ocean  basins  as  well. 
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The  forecasting  scheme  presented  in  this  study  may  not  significantly 
exceed  the  accuracy  of  the  forecasts  produced  by  the  experts  at  the 
National  Hurricane  Center  or  the  Joint  Typhoon  Warning  Center  but  it 
does  give  quantitative  guidance  to  the  parameters  which  appear  to  be 
important  to  intensity  change.  Because  the  scheme  was  designed  for  use 
on  standard  isobaric  charts  using  available  information,  with  minimal 
training  any  forecaster  should  be  able  to  make  the  computations  in  any 
weather  station  in  only  10-15  minutes.  Also,  the  selection  of  the  850 
and  200  mb  charts  allows  us  to  supplement  the  wind  data  with  satellite 
derived  winds  to  increase  our  data  coverage.  The  scheme  can  also  be 
programmed  for  computer  printout. 

Although  this  study  failed  to  predict  the  correct  intensity  change 
in  approximately  25%  of  the  cases,  the  scheme  performs  as  well  as  cur¬ 
rent  forecasting  techniques.  You  will  always  have  an  Eloise  to  defy 
your  best  physical  reasoning.  However,  the  idea  that  intensity  change 
can  be  approached  in  a  quantitative  manner  bears  further  investigation. 
More  research  should  be  done  to  determine  how  the  physical  coupling 
between  tropical  storms  and  large  scale  baroclinic  zones  takes  place. 
Hopefully,  some  of  the  ideas  presented  in  this  study  will  stimulate 
further  research  in  an  attempt  to  better  understand  intensity  change 
processes. 
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APPENDIX  A 


Appendix  A  is  a  summary  of  the  parameters  analyzed  for  all  storms 
which  make  up  this  study.  Tables  7-10  represent  the  data  for  the  orig¬ 
inal  composite  study.  Tables  11-14  list  all  results  for  the  independent 
storms  which  were  tested.  All  parameters  analyzed  (except  6°  belt  ^qq) 
represent  the  5-11°  belt  average  over  octants  8,  1  and  2.  Units  for  the 
respective  parameters  are  as  follows: 


zonal  wind  parameters  (Ug50>  U500»  u2oo^ 
vertical  zonal  wind  shear  (U85()-U500,  U^-U^) 
temperature  field  (Tg^) 
height  field  (H^) 
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DEP  represents  a  storm  in  the  depression  stage  (max  wind  <_  34  kts) .  Prior  to  1974  the  depressions  were 
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APPENDIX  B 

COMPUTATIONAL  PROCEDURES  FOR  THE  INTENSITY  CHANGE  FORECASTING  PARAMETER 


The  computation  of  the  intensity  change  parameter  (IP)  has  been 
designed  so  that  it  could  be  performed  quickly  and  with  data  available 
at  any  weather  station.  It  uses  the  standard  isobaric  850,  500,  and 
200  mb  charts  which  are  available  twice  daily.  With  only  a  minimum  of 
practice  the  IP  could  be  computed  in  less  than  15  minutes  using  the 
following  guidelines: 

1.  Locate  the  storm  center  on  the  850,  500  and  200  mb  charts 
based  on  the  position  given  by  the  Miami  Hurricane  bulletin 
(or  the  bulletin  issued  by  JTWC)  for  a  particular  time  period. 

2.  Using  a  compass  mark  off  the  6,  8,  and  10°  radial  belts  on  the 
north  side  of  the  storm. 

3.  Subdivide  these  radial  belts  into  the  45°  octants  8,  1  and  2 
as  in  Fig.  2.  For  storms  affecting  the  eastern  U.S.,  China, 
Korea  or  Japan  use  octants  1,  2  and  3. 

4.  Compute  average  values  of  the  parameters  Tg5Q,  H2oo>  u850» 

U500  and  U200  over  the  5-11°  belt  over  octants  8,  1  and  2 
(or  1,  2  and  3)  by: 

a.  Obtain  an  average  value  in  each  octant  for  each  of  the 
5-7°,  7-9°  and  9-11°  radial  belts. 

b.  Average  the  values  from  each  octant  in  each  radial  belt  to 
get  three  belt  averages. 

c.  Average  these  three  belt  averages  together  and  the  resul¬ 
tant  number  represents  the  5-11  belt  average  over  the 
desired  octant. 

A  quick  streamline  analysis  may  have  to  be  done  to  make  it  easier  to 
resolve  the  zonal  wind  components.  When  the  wind  components  have  been 
determined  in  knots,  divide  by  2  to  get  the  wind  speed  in  1  s  \  It  is 
helpful  to  use  a  table  such  as  those  found  in  Appendix  A  to  keep  track 
of  all  parameters. 

5.  Calculate  the  850  to  500  mb  and  850  to  200  mb  wind  shears  from 
the  table. 
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6.  Using  either  the  graph  in  Fig.  44  or  45  plot  each  of  the  final 
parameters  on  the  appropriate  scale  and  read  off  the  value  from 
the  bottom  of  the  graph.  Add  these  five  numbers  together  and 
multiply  by  10.  The  resultant  number  is  the  intensity  change 
parameter  (IP). 

Further  refinements  can  be  made  in  the  Atlantic  basin  based  on 
the  seasonal  variations  described  in  section  4.2.  The  adjustments  to  be 
made  to  IP  are  given  in  Table  15.  After  all  computations  and  adjustments 
have  been  made,  if  IP  >  0  we  should  expect  an  intensifying  system. 
Negative  values  would  represent  a  filling  or  steady  storm. 


TABLE  15 

Correction  factors  to  be  applied  to  IP  based  on  seasonal  variation. 
Apply  the  appropriate  correction  based  on  date  of  occurrence  of  the 
cyclone  system. 


TIME  PERIOD 

MAY 

15-31 

JUN 

1-15 

m3 

H2 

JUL 

1-15 

JUL 

16-31 

AUG 

1-15 

AUG 

16-31 

SEP 

1-15 

[^2 

OCT 

1-15 

OCT 

15-31 

ADJUSTMENT 

FACTOR 

+80 

+  50 

+20 

-10 

-30 

-50 

-50 

-30 

-10 

+20 

+50 

Computation  graph  for  IP  in  the  Atlantic  basin 


Computation  graph  for  IP  in  the  Pacific  basin. 
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APPENDIX  C 

DESCRIPTION  OF  COMPOSITE  DATA  SETS 

Dropco  Intensifying;  Tropical  storms  and  hurricanes  which  either  form¬ 
ed  or  moved  into  the  Gulf  of  Mexico  and  showed  significant  intensifica¬ 
tion  (wind  speed  change  >20  kts  in  24  hours)  prior  to  landfall  along 
the  U.S.  gulf  coast.  The  data  came  from  the  official  best  track  posi¬ 
tion  of  the  National  Hurricane  Center  for  the  years  1957-1977  and  in¬ 
clude  22  cases.  The  stratification  was  based  primarily  on  wind  speed 
change  and  central  pressure  when  available.  The  composite  was  run  for 
two  time  periods:  42  hours  prior  to  landfall  with  average  latitude  = 
23.3°N,  longitude  =  88.4°W  and  18  hours  prior  to  landfall  with  latitude 
*  25.1°N,  longitude  =  90.7°W. 

Dropco  Non- intensifying:  Tropical  storms  and  hurricanes  which  either 
formed  or  moved  into  the  Gulf  of  Mexico  and  showed  little  or  no  ten¬ 
dency  for  intensification  prior  to  landfall.  The  data  were  selected 
from  the  best  track  position  of  the  National  Hurricane  Center  for  the 
years  1957-1977  and  include  13  cases.  The  stratification  was  made  as 
in  the  intensifying  case  and  was  also  run  for  two  time  periods:  42 
hours  prior  to  landfall  with  latitude  =  23.8°N,  longitude  =  84.7°W  and 
18  hours  prior  to  landfall  with  latitude  =  25.4°N,  longitude  =  88.0°W. 

Holliday  (1979)  rapidly  Deepening:  Storms  from  the  Northwest  Pacific 
for  the  years  1961-1970  where  24  hour  deepening  rate  was  42  mb  day  . 
The  data  were  stratified  into  three  periods:  24  hours  prior  to  the  onset 
of  rapid  deepening,  the  period  at  onset,  and  24  hours  after  the  onset 


of  deepening. 


1 


Bill  Frank  (1976)  Pacific  Data  Sets:  These  data  sets  were  made  from 
ten  years  (1961-1970)  of  Northwest  Pacific  rawinsonde  data.  The  data 
were  stratified  by  both  latitude  and  tendency  in  the  following  manner: 


Type  I 


Type  II 


P  <  980  mb 
c  — 


980  mb  <  P  <  1000  mb 
c  — 


Type  III  P^  >  1000  mb 


Type  A 


Storm  which  at  some  point  reaches  P^  <^  970  mb 


Type  A  Early  Same  as  Type  A  but  980  mb  <  Pc  £  1000  mb  and  in  a 
deepening  trend  towards  a  typhoon. 


Type  B 


Storm  with  980  mb  <  P  <_  1000  mb  but  P  never  drops 
below  980  mb. 


lp  represents  central  pressure  of  the  storm, 
c 


Bill  Frank  (1976)  West  Indies  Data  Sets:  These  data  sets  were  made 
from  14  years  (1961-1974)  of  western  Atlantic  rawinsonde  data.  The 
data  were  stratified  by  latitude  and  filling  or  deepening  tendencies  as 
taken  from  the  best  track  positions  of  the  National  Hurricane  Center. 


Zehr  (1976)  Developing  Trade  Wind  Clusters:  Pre-typhoon  disturbances 
found  in  the  trade  wind  belt  taken  from  the  time  they  were  first  identi¬ 
fied  as  cloud  clusters  until  they  attained  an  intensity  with  maximum 
sustained  winds  of  50  knots.  The  data  were  taken  from  the  central  and 
western  North  Pacific  during  the  ten  year  period  1961-1970. 
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Zehr  (1976)  Stage  0  Non-developing  Cluster:  Non-developing  cloud  clus¬ 
ters  in  the  central  and  western  Pacific  for  the  period  1967-1968.  The 
clusters  were  found  in  the  region  from  the  equator  to  18°N  and  125°E  to 
160°W. 


